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Abstract 

We investigate the Wightman function, the vacuum expectation values of the field squared 
and the energy-momentum tensor for a massless scalar field with general curvature coupling 
parameter in spatially flat Friedmann-Robertson- Walker universes with an arbitrary number 
of toroidally compactified dimensions. The topological parts in the expectation values are 
explicitly extracted and in this way the renormalization is reduced to that for the model 
with trivial topology. In the limit when the comoving lengths of the compact dimensions 
are very short compared to the Hubble length, the topological parts coincide with those 
for a conformal coupling and they are related to the corresponding quantities in the flat 
spacetime by standard conformal transformation. This limit corresponds to the adiabatic 
approximation. In the opposite limit of large comoving lengths of the compact dimensions, 
in dependence of the curvature coupling parameter, two regimes are realized with monotonic 
or oscillatory behavior of the vacuum expectation values. In the monotonic regime and 
for nonconformally and nonminimally coupled fields the vacuum stresses are isotropic and 
the equation of state for the topological parts in the energy density and pressures is of 
barotropic type. For conformal and minimal couplings the leading terms in the corresponding 
asymptotic expansions vanish and the vacuum stresses, in general, are anisotropic though 
the equation of state remains of barotropic type. In the oscillatory regime, the amplitude 
of the oscillations for the topological part in the expectation value of the field squared can 
be either decreasing or increasing with time, whereas for the energy-momentum tensor the 
oscillations are damping. The limits of validity of the adiabatic approximation are discussed. 



1 Introduction 

Friedmann-Robertson- Walker (FRW) spacetimes are among the most popular backgrounds in 
gravitational physics. There are several reasons for this. First of all current observations give a 
strong motivation for the adoption of the cosmological principle stating that at large scales the 
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Universe is homogeneous and isotropic and, hence, its large-scale structure is well described by 
the FRW metric. Another motivation for the interest to FRW geometries is related to the high 
symmetry of these backgrounds. Due this symmetry numerous physical problems are exactly 
solvable and a better understanding of physical effects in FRW models could serve as a handle 
to deal with more complicated geometries. In particular, the investigation of quantum effects 
at the early stages of the cosmological expansion has been a subject of study in many research 
papers (see, for example, [HE] and references therein). The original motivations in this research 
were mainly to avoid the initial singularity and to solve the initial conditions problem for the 
dynamics of the Universe. The further increase of the interest to this topic was related to the 
appearance of inflationary cosmological scenarios |3j. During an inflationary epoch, quantum 
fluctuations introduce inhomogeneities and may affect the transition toward the true vacuum. 
These fluctuations provide a mechanism which explains the origin of the primordial density 
perturbation needed to explain the formation of the large-scale structure in the universe. 

In many problems we need to consider the physical model on the background of manifolds 
with compact spatial dimensions. In particular, the idea of compact extra dimensions has been 
extensively used in supergravity and superstring theories. From an inflationary point of view, 
universes with compact dimensions, under certain conditions, should be considered a rule rather 
than an exception The models of a compact universe with non-trivial topology may play 
important roles by providing proper initial conditions for inflation. The quantum creation of 
the universe having toroidal spatial topology is discussed in [5] and in references (Sj within the 
framework of various supergravity theories. In the case of non-trivial topology, the boundary 
conditions imposed on quantum fields give rise to the modification of the spectrum for vacuum 
fluctuations and, as a result, to the Casimir-type contributions in the vacuum expectation values 
of physical observables (for the topological Casimir effect and its role in cosmology see [8] and 
references therein). In the Kaluza-Klein-type models, the topological Casimir effect has been 
used as a stabilization mechanism for moduli fields which parametrize the size and the shape of 
the extra dimensions. Recently, the relevance of the Casimir energy as a model for dark energy 
needed for the explanation of the present accelerated expansion of the universe has been pointed 
out (see [9] and references therein). 

In most work on the topological Casimir effect in cosmological backgrounds, the results for 
the corresponding static counterparts were used replacing the static length scales by comoving 
lengths in the cosmological bulk. This procedure is valid in conformally invariant situations or 
under the assumption of a quasi-adiabatic approximation. For non-conformal fields the calcula- 
tions should be done directly within the framework of quantum field theory on time-dependent 
backgrounds. In a series of recent papers we have considered the topological Casimir effect for 
a massive scalar field with general curvature coupling [TO] and for massive fermionic field [H] 
in background of the de Sitter spacetime with toroidally compactified spatial dimensions. Due 
to the high symmetry of the de Sitter spacetime in both cases the corresponding problems are 
exactly soluble. 

In the present paper we will consider another exactly soluble problem for the topological 
Casimir effect on background of a more general class of spacetimes. Namely, we investigate the 
vacuum expectation values of the field squared and the energy-momentum tensor for a massless 
scalar field with arbitrary curvature coupling parameter induced by the compactness of spatial 
dimensions in spatially flat FRW universes where the scale factor is a power of the comoving time. 
The vacuum polarization and the particle creation in the FRW cosmological models with trivial 
topology have been considered in a large number of papers (see [H HI [7] and references therein 
for early research and [12j-|23j and references therein for later developments). In particular, the 
vacuum expectation values of the field squared and the energy-momentum tensor in models with 
power law scale factors have been discussed in [UJ [15], [M]- [22] (see also [2]). The expectation 
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value of the field squared is the crucial quantity in discussions of the topological mass generation 
as a result of one-loop quantum corrections and dynamical symmetry breaking. It also is of 
interest in considerations of phase transitions in the early universe, in particular, related to 
inflationary scenarios. In addition to describing the physical structure of the quantum field at a 
given point, the expectation value of the energy-momentum tensor acts as the source of gravity 
in the Einstein equations. It therefore plays an important role in modelling a self-consistent 
dynamics involving the gravitational field. 

The present paper is organized as follows. In the next section we evaluate the positive- 
frequency Wightman function in spatially flat FRW model with topology BP x [S 1 ) q . By using 
the Abel-Plana summation formula, we decompose this function into two parts: the first one 
is the corresponding function in the geometry of the uncompactified FRW spacetime and the 
second one is induced by the compactness of the spatial dimensions. In section [3] we use the 
Wightman function for the evaluation of the vacuum expectation value of the field squared. The 
asymptotic behavior of the topological part is investigated in the early and late stages of the 
cosmological evolution. In section [5] we consider the vacuum expectation value of the energy- 
momentum tensor. The part in this expectation value corresponding to the uncompactified FRW 
model is well investigated in the literature and we are mainly concerned with the topological 
part. The special case of the model with a single compact dimension and the corresponding 
numerical results are discussed in section [5l The main results of the paper are summarized in 
Section [6j 



2 Wightman function 

We consider a scalar field with curvature coupling parameter £ evolving on background of the 
(D + l)-dimensional spatially flat FRW spacetime with power law scale factor. Such a field is 
described by the equation 

(ViV 1 + m 2 + ^ = 0, (1) 

where V; denotes the covariant derivative and R is the curvature scalar of the background 
spacetime. For the most important special cases of minimally and conformally coupled fields 
one has the values of the curvature coupling parameter £ = and £ = £d = (D — 1)/4D 
respectively. The spatially flat FRW line element expressed in comoving time coordinate is 

D 

ds 2 = dt 2 - a 2 (t) ^(d^) 2 , a(t) = at c , (2) 
i=i 

with the curvature scalar 

R = Dc[(D + l)c-2]/t 2 . (3) 

From the (D + l)-dimensional Einstein equations we find the energy density p and the pressure 
p for the source of the metric corresponding to ([2]): 

_ D{D - 1) c 2 _ D - 1 c(2 - cD) 
P ~ 16vrG P ,P ~16^G J2 > ( 4 ) 

where G is the (D + l)-dimensional gravitational constant. The corresponding equation of state 
has the form p = [2/(cD) — l]p. In the special case 

c = 2/(D + l), (5) 

the equation of state is of the radiation type and the Ricci scalar vanishes. The case c = 2/D cor- 
responds to the dust-matter driven models. Note that the line element ([2]) with c > 1 describes 
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the power law [33] and extended [31] inflationary models. Such a power law expansion may 
be realized if a scalar field <fi with an exponential potential V((f>) = dominates the energy 
density of the universe. The corresponding cosmological models admit power law inflationary 
solutions of the form ([2]) (see, for instance, [35] ) with c = 4A -2 / (D — 1). Solutions having power 
law scale factor with c = 2(n — l)(2n — 1)/(D — 2n) arise in higher-order gravity theories with 
the Lagrangian R n [36] . 

For the further discussion, in addition to the synchronous time coordinate t it is convenient 
to introduce the conformal time r in accordance with dt = a(t)dr or 

t=[a(l-c) T ] 1 ^ 1 - c \ (6) 

Here we assume that c / 1. In the case c = 1 one has t = toe aT , with a constant to, and this 
case should be considered separately. In terms of the conformal time the line element takes the 
form: 

D 

ds 2 = n 2 (r)[dT 2 - ^(dz 1 ) 2 }, «(r) = a[a(l - c)t\ c ^ x - c \ (7) 
i=l 

which is manifestly conformal to the Minkowski spacetime. Note that one has ^ r < oo for 
< c < 1 and — oo < r ^ for c > 1. 

We will assume that the spatial coordinates z l , I = p + 1, . . . ,D, are compactified to S 1 : 
^ z l ^ Li, and for the other coordinates we have — oo ^ z l ^ +oo, / = 1, . . . ,p. Hence, we 
consider the spatial topology BP x (S 1 ) q with p + q = D. For p = 0, as a special case we obtain 
the toroidally compactified FRW spacetime. Along the compact dimensions we will consider the 
boundary conditions 

<p(r, z p ,z q + e t Li) = e 2mai tp(r, z p , z q ), (8) 

where z p = (z 1 , . . . , z p ), z q = (z p+1 , . . . , z D ), and e^, I = p + 1, . . . , D, is the unit vector along 
the direction z . In ([8]) we have introduced phases with constants ^ a\ ^ 1. The special 
cases a\ = and a\ = 1/2 correspond to periodic and antiperiodic boundary conditions. The 
corresponding fields are also called as untwisted and twisted ones. 

Among the most important characteristics of the vacuum state are the vacuum expectation 
values (VEVs) for the field squared and the energy-momentum tensor. These VEVs are obtained 
from the Wightman function and its derivatives in the coincidence limit of the arguments. For 
the topology under consideration we will denote this function by Gp q (x,x'). The Wightman 
function is also important in the consideration of the response of particle detectors (see, for 
instance, [I]). Let {(p a (x), ip* a (x)} be a complete set of solutions to the classical field equation 
satisfying the periodicity conditions ([8]) and a is a set of quantum numbers specifying the 
solution. The Wightman function may be evaluated by using the mode-sum formula 

G+ q (x,x') = (Ol^MxOlO) = J> CT (xK(x'), (9) 

where |0) is the amplitude for the vacuum state under consideration. 

In accordance with the symmetry of the problem the mode solutions of the field equation are 
separable: (f a (x) = 0(r)e ikp Zj ' +lk ' J ' Z9 , where k p = {k\, k 2 , . . . , k p ) and k q = (k p+ i,k p+ 2, . . . , &d). 
For the components of the wave vector along the uncompactified directions we have — oo < 
k\ < oo, I = 1,2, ... ,p. The components along the compact dimensions are quantized by the 
periodicity conditions ([8|) and for the corresponding eigenvalues we have 

k = 2vr(n z + a{)/L u n t = 0, ±1, ±2, . . . , I = p + 1, . . . , D. (10) 
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From the field equation ([T]) we obtain an ordinary differential equation for the function 0(r): 



** + + + a' { [„•(! - cM */<-W + { PC ,"° +1 >7 21 ^ = 0, (11) 

ar z (1 — cjr ar [ [a(l — c)t\ z 



fc=,/k2+ k 2 (12) 



where 



For a massless scalar field, m = 0, the general solution of this equation is expressed in terms of 
a combination of Hankel functions: 



c 1 Hi 1 \kT)+c 2 HP(kT) , (13) 



with the notations 



and 

1 



^ = ^ h = W=Ty (14) 



i V(cO-l) 2 -^c[(D + l)c-2]. (15) 

Note that ^ is either real and nonnegative or pure imaginary. For a conformally coupled field 
we have v = 1/2 and for a minimally coupled field v = |6|. Different choices of the coefficients 
in (|13|) correspond to different vacuum states. We will consider the Bunch-Davies vacuum [26] 
for which c\ = 0. With this choice, in the limit c — > we recover the standard vacuum in the 
Minkowski spacetime associated with the modes k~ 1 / 2 exp(ik p • z p + ik. q ■ z q — ikr). 
The corresponding eigenfunctions satisfying the periodicity conditions take the form 

<p a (x) = C^H^ikT)^^"- 2 ". (16) 

The coefficient C a with a = (k p ,n p +i, . . . , tid) is found from the orthonormalization condition 

d D xy/\g\g 00 [ip a {x)d T ip* a ,{x) - tp*,(x)d T (p a (x)] = i5 a(7 >, (17) 



where the integration goes over the spatial hypersurface r = const, and 8 aa i is understood as 
the Kronecker delta for discrete indices and as the Dirac delta-function for continuous ones. By 
using the well known properties of the Hankel functions, this leads to the result 

where and in what follows we denote by 

V q = L p+1 ■■■L D (19) 

the volume of the compact subspace. 

In the case p = 0, q = D (topology (S 1 ) ) and ai = 0, I = 1, . . . , D, there is also a normaliz- 
able zero mode with kp = which is independent of the spatial coordinates (for the discussion 
of zero modes in topologically nontrivial spaces see [37] ). The corresponding eigenfunction has 
the form </?o( T ) = V b ( c ioV l/ + c 2O r 7~ I 0- From the normalization condition we have the relation 
c*qC2o — C10C20 = Mgn(r)/(2az/Vb) for real v and the relation | C20 1 2 — |cio| 2 = sign(r)/(2|i/|aVb) 
for imaginary v. These formulae give a relation between the coefficients c\q and C20 and one of 
them remains undetermined. This arbitrariness arises because of freedom to choose the quan- 
tum state |37j . In the discussion below we will consider the contribution to the VEVs from the 
nonzero modes. If the zero mode is present, its contribution should be added separately. 
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Substituting the eigenfunctions (|16p with the normalization coefficient (|18|) into the mode- 
sum formula for the Wightman function, one finds 



g±jx,x') = A{j1 i)l ! dk p j k »- Az * y 

p ' qV ' 2PirP +1 V n ^ 



x K v {k V e signiT)ni/2 )K„ (kr]'e- s[gn{T)ni/2 ), (20) 
where n q = (n p+ i, . . . , ud), Az l = z l — z l> , and 

A = a 1 - D [a\l-c\]( D - 1 ^ c ~ 1 \ (21) 



In formula (|20p . for the further convenience, we have written the Hankel function in terms of the 
modified Bessel function. For a scalar field with periodic boundary conditions along compact 
dimensions (ai = 0) and for Hev ^ p/2 the integral in (|20p has an infrared diverg ence at k p — 
coming from the modes n q = 0. In this case the Bunch-Davies vacuum is well defined only if 
Rev < p/2. For the case p = and a\ = 0, in ([20|) n q = rip ^ and the contribution of the 
additional zero mode should be added. 

We apply to the sum over n p+ \ in (|20p the Abel-Plana type summation formula 

V g(n + /3)f(\n + /3\) = / du [g{u) + g{-u)] f (u) 

n=-oo J ° 

+* r du [f(iu) - K-iu)] y ^S^hi - (22) 

1/0 A=±l 

This formula is obtained by combining the summation formulae given in |38j . After the ap- 
plication of this formula, the term in the expression of the Wightman function with the first 
integral on the right of (|22p corresponds to the Wightman function in the FRW model with p + 1 
uncompactified and q — 1 toroidally compactified dimensions, G^ +1 q _ 1 (x,x'). As a result one 
finds the recurrence formula 

G+ 9 (x, x') = G+^.^x, x') + A p+1 G+ q (x, x'), (23) 

where the second term on the right is induced by the compactness of the z p+1 - direction and is 
given by the expression 



A ^G+ (r r') - A (w') b f n ik p -Az p 

Ap+l^ p>q {X,X ) - 2(27r)p+ l^_ i J dk P e 2^ 



e ik 9 _i-Azq_i 



dy - I y y 

Jy 2 + k 2 , + k^ =i t^x e Lp + 1 v /2/2+k ^ +k "«-i +27raQp+1 - 1 
x { ff„ ( w ) [/_„ (7/y) + 4 (r/'y )] + [I„ (m) +I-u(vy)} K v ( V 'y) } , (24) 
where I u (z) is the modified Bessel function. Here the notations V^_i = £ p +2> • • • , and 

D 

n q -i = (n p+2 , ■ ■ ■ , n D ), k 2 nq l = y (2irni/Li) 2 , (25) 

l=p+2 
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are introduced. In formula (|24p the integration with respect to the angular part of k„ can be 
done explicitly. By using the recurrence formula (|23p . the Wightman function is presented in 
the form 

D-l 

G+ q (x,x') = G+ RW (x,x') + £ A i+1 G+ D _.(x,x'), (26) 

j=p 

where Gp RW (x,rc') = Q {x,x') is the corresponding function in the FRW model with spatial 
topology R D and the second term on the right is the topological part. 

The integral representation for the function Gp RW (x, x') is given by Eq. ([20]) with p = D 
and q = 0. The corresponding integral is explicitly evaluated by using the formula from |39| and 
we find 

G+ (x xl) _ A(w')^ T{D/2 + v )T(D/2 - u) p (i-D)/2 ( ) m 

where P^{u) is the associated Legendre function of the first kind and 

n= (|Az| 2 -r ? 2 -?? /2 )/(2?7?7 / ). (28) 

An alternative representation is obtained by making use of the relation between the associated 
Legendre function and the hyper geometric function (see, for instance, [40j ) : 

A(W)^T(D/2 + unD/2-u) (D D D + I 1 - u\ 
^ FRW [x,x)- (47r)(D+1)/2 r((£) + l)/2) V 2 + ^' 2 2 ' 2 J ' {ZJ) 

Note that we have 

A( m ')^ = \1- c\ D -\l£)^- D ^ 2 . (30) 

In the special case D = 3 and £ = the formula (|29p coincides with the result given in [26J (with 
the missprint in the coefficient corrected, see also PQ). Shifting the time coordinate t — » t + cto, 
with a constant io an d taking the limit c — > oo for fixed t and to, from ([2]) we obtain the 
corresponding line element in de Sitter spacetime with the scale factor a(t) = e 4// '°. In this limit 
the expression (jUD tends to t^ D and from (EHJ) we recover the expression for the Wightman 
function in (D + l)-dimensional de Sitter spacetime. 



3 VEV of the field squared 

The VEV of the field squared in FRW model with topology R p x (S 1 ) 9 is obtained from the 
two-point function GfJx,x') taking the coincidence limit of the arguments. In this limit the 
Wightman function diverges and a renormalization procedure is necessary. The crucial is that 
under the toroidal compactification the local geometry is not changed and the divergences are 
contained in the part corresponding to the uncompactified FRW spacetime. As we have already 
extracted this part, the renormalization is reduced to that in the uncompactified FRW model 
which is well investigated in the literature (see for example [U [2] and references given therein) . 
For the VEV of the field squared we have the recurrence formula 

(^ 2 )p,<? = (y 2 ) P +i l(7 -i + A p+ i(</? 2 ) Pi(? , (31) 

where (</? 2 ) p +i, g -i is the expectation value in the FRW spacetime with the spatial topology 
R p+1 x (S' 1 ) 9-1 and the second term on the right is the part due to the compactness of the 
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z p+1 - direction. The latter is directly obtained from (|24p taking the coincidence limit of the 
arguments. The integral over k p is evaluated by making use of the formula 

( k p + a 2 y 1/2 _ A icy„\( p -i)/2 ^ cos(27rna p 

/ ^ „-Lti-H 

A=±l ' 

with the notation f v (x) = x v K v {x). Now for the topological part in (|31H we find 



£ r ^ sggr — = 4 ( 2.,^')/2 E ~ym^o ( a) (32) 



AArj 2b \ f°° 

(2vr)^V Vi n! ^ Zrl io 



This integral representation is valid for Rei^ < 1. 

Similar to (|26p . we find the following decomposition for the VEV of the field squared: 

D-l 

( V 2 ) m = <v? 2 )frw + (<f?)®, {V 2 )t\ = E A i+i<<A^-i> (34) 

where (v? 2 )frw = (^ 2 ) D,o is the VEV in the spatial topology R D and the part (<p 2 ) P % is induced 
by the nontrivial topology. As it can be seen from formula ([33]), the topological part has the 
following general structure: {f 2 )p^q = t 1- ® f(L p+ i/rj, . . . ,Ld/t]), where the form of the function 
f(x±, . . . , x q ) is directly obtained from (|33|) and (|34|) . The ratios in the arguments of this 
function may also be written as 

L l /r ] = {\l-c\/c)Lf ) /r m (35) 

with = QLi and rn = t/c being the comoving length of the compact dimension and the 
Hubble length respectively. 

For a conformally coupled massless scalar field one has v = 1/2 and [I- U (x) + I u (x)] K v [x) = 
1/x. The integral in (|33p is explicitly evaluated and we find 

a / 2\ 2Arj 2b - 1 sr^ ^ cos(2vrna p+ i) 

A p+ i(y ) p , g = w ^ E ^7~^ f P /2(nL p+1 k nq ^). (36) 



P+l 



We could obtain this result directly from the corresponding formula in the Minkowski spacetime 
with topology R p x (S' 1 )' 2 , by taking into account that Ar] 2b ~ l = Q 1 ~ D and by using the fact that 
two problems are conformally related: {(p 2 )p\ = {(p 2 )p t 'q . The expression for {f 2 ) p \' q M ^ is 
directly obtained from (|36p and is valid for arbitrary values of the curvature coupling parameter 
£. Note that we have v = 1/2 in the general case of the curvature coupling parameter for 
the power law expansion with c = 2/{D + 1) which corresponds to radiation driven models. 
In this case the topological part of the VEV is given by formula (|36p with b = —1/2 and 
A = a- D ~ 1 {D + 1) 2 /{D-1) 2 . 

Now we turn to the investigation of the topological part A p+ i(<^ 2 ) Pj(? in the VEV of the 
field squared in the asymptotic regions of the ratio L p+ i/r]. For small values of this ratio, 
Lp+i/f] *C 1, we introduce the new integration variable y = L p+ \x. By taking into account 
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1/x, we find that to the leading order 
An+i(9? 2 )p i(? coincides with the corresponding result for a conformally coupled massless field: 



that for large values x one has [I- U (x) + I v {x)\ K v (x) 



2n 



l-D 



(2vr)P/ 2 + 1 y 9 _i 



E E 

_ 1 £Zi~ 1 n=l 



cos(2-7rna p+ i) 
(nL p 



f p j 2 (nL p+1 k nq -i), L p+1 /rj < 1. (37) 



In terms of the comoving time coordinate the topological part behaves as A p+ i{(p 2 ) Ptq oc t^^ D ^ c . 
The limit under consideration corresponds to early stages of the cosmological expansion (i —* 0) 
for c > 1 and to late stages (t — ► oo) for < c < 1. Note that, by taking into account 
([35]) . the condition L p+ \/rj <C 1 can also be written in the form <C rn- Hence, the 

asymptotic formula (|37p corresponds to small proper lengths of compact dimensions compared 
with the Hubble length. The expression on the right hand side of (137ft is obtained from the 
corresponding expression in Minkowski spacetime with the topology R p x {S l ) q and with the 

(c) 

lengths of the compact dimensions L p +\, . . . ,Ld, replacing L\ by the comoving lengths . 
Hence, (|37p corresponds to the adiabatic approximation. In this regime the topological part 
monotonically decreases with increasing scale factor. 

In the opposite limit of large values for the ratio L p+ i/r], the behavior of the topological part 
is qualitatively different for real and imaginary values of the parameter v. For positive values 
v, by using the asymptotic formulae for the modified Bessel functions for small values of the 
arguments, to the leading order we find 



v+l A„2b-2v 



l Al] 



(2tt)(p+ 3 )/ 2 Vi 



E E 



cos(27rna p+ i) 

~2v 



Xf(p+l)/2-A nL P+l k n q -l), L p+1 /r] > 1. 



(38) 



In terms of the synchronous time coordinate we have A p+ i(ip 2 ) Ptq oc t 2 ( c - 1 ) u - cD + 1 _ This limit 
corresponds to early stages of the cosmological expansion (t — > 0) for < c < 1 and to late 
stages of the cosmological expansion (t — > oo) when c > 1. Note that for < c < 1/D and £ < 
one has b > u and the topological part in the expectation value of if 2 increases with increasing 
comoving time t. For a minimally coupled scalar field and for c G (0, l/D)Li(l, oo) one has v = b, 
and in the limit under consideration the VEV tends to finite nonzero value. In the case v = 
the corresponding asymptotic formula is obtained from (|38|) replacing T(v) — ► 41n(L p+ i/r/) and 
after substituting v = 0. The limit under consideration corresponds to large comoving length of 
the compact dimension compared to the Hubble length: L p \_-y 3> rn- 

For pure imaginary values v and in the same limit L p+ \/r] S> 1, we use the formula 



K u (z) [I_ v ( z ) + /„(*)] ~ Re 



rfr)(*/2) 
rn - 1/ 



-2u 



valid for small values z. Now we can see that in the leading order 

a / 2v ^AATj 2b c ^[2\v\\n{L p+1 /r]) + <j>] 

A p+1 (^> ) M « (27r) ( P+3 )/2 L P +iyg ■ L p+i/^ » 

where the constants B and <j> are defined by the relation 



(39) 



(40) 



n=l 



(41) 



n„_ieZ9-l 
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Hence, the behaviour of the topological part in the VEV for the field squared is oscillatory in the 
early (late) stages of the cosmological expansion for < c < 1 (c > 1). The amplitude of these 
oscillations behaves like t 1-cD and the distance between the nearest zeros linearly increases with 
increasing t. Note that for < c < 1/D the oscillation amplitude increases with time. 

From the asymptotic analysis given above it follows that in the special case of a minimally 
coupled field the topological part A p+ i((p 2 ) Pjq tends to finite nonzero value at early stages of 
the cosmological expansion when < c < 1/D. At late (early) stages and for c < 1 (c > 1) the 
topological part behaves like t^~ D ^ c . At early times and for 1/D < c < 1 one has A p+ i(ip 2 ) p q oc 

t 2(l-cD) > 



4 VEV of the energy- momentum tensor 

Having the Wightman function and the VEV of the field squared, the expectation value of the 
energy-momentum tensor is evaluated with the help of the formula 



(T ik ) m = Jim did' k G+ q {x,x') + U--^j feV ; V ; - £V,V fc - £R ik 



(V 2 ) p , q , (42) 



where Ri k is the Ricci tensor for the spatially flat FRW spacetime with the components 

Dct~ 2 c(Dc-l) 
Roo = 7-1 Rik = —$ikTi vT~ 9 ' i,k = 1, . . . ,D. (43) 

C — 1 (1 — CPT Z 

In Eq. (|42p we have used the expression of the metric energy-momentum tensor for a scalar 
field which differs from the standard one (see, for example, [I]) by the term vanishing on the 
solutions of the field equation |41j . The topological part in the VEV of the energy-momentum 
tensor is obtained substituting (I23p and (|3ip into formula (142p . In this calculation we need the 
covariant d'Alembertian acted on the topological part of the field squared: 

^^i/A , 2y 2 [1 -^l 2 Ac ^ ^ cos(2vrna p+ i) f°° , 3 _ 2b 

ng_i<EZi y n— x 

xf( P -i)/2(nL p+1 ^y 2 + k2 ? _ x ) I u (z)K u (z)\ z=m . (44) 

Here and in the discussion below, to simplify our notations, we have introduced the functions 

K v {z) = z b K p (z), I v (z) = z b [I p (z) + /_„(*)] , (45) 

and the function f u (x) is defined after formula (f32j) . 

For the VEV of the energy-momentum tensor we have the following recurrence formula 

(T t % q = (Tf) P+ i, q -i + A p+1 (T/%, 9 , (46) 

where the second term on the right hand side is due to the compactness of the direction z p+1 . 
For the topological part in the energy density we find 

A /T°\ Am ~ 2 cos(27ma p+1 ) f°° 3 -2fe 

p+i( o)p ' q ~ p^^i^fI dyy 

x f( P -m(nL p+1 ^y 2 + k^)F (0) (w), ( 47 ) 

n 9 _iSZ9-i 
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with the notation 



F®(z) = \l' u {z)K> u {z) + -£&(I v K v y - \ 
I z{\ — c) 2 



£D{D - 1)6 
z 2 (l- 



.21 



I V K V . 



(48) 



In a similar way, for the vacuum stresses the following representation holds (no summation over 





\+l( T i)p,< 



4AQ- 



( 27r )( P +3)/2^ 



9-1 



E E 



cos(27rna p+ i) 



dyy 



3-26 



(49) 



2=nL p+ W?/ 2 +k2 (; _ 1 



where we have used the notations 



4 p+1) (y) 



f(p+i)/2(y), i = l,2,...p, 

-pf( P +i)/2(y) - y 2 f( P -i)/2{y), 

, p \»j - (nL p+ ik i ff (s _ 1)/2 (y),i=p+2,...D. 
In formula (1491) . the function F(z) is defined as 



(50) 



F(z) 



+2 



1 _ (Dc-2)(g-g D )+^ 
? 4 ? z 2 (l-c) 2 



I v {z)K v {z) 



(51) 



Now it can be checked that the topological parts in the VEVs of the field squared and the 
energy-momentum tensor obey the trace relation 



(52) 



In particular, for a conformally coupled field the corresponding energy- momentum tensor is 
traceless. 

The recurrence relation (|46p allows to present the VEV of the energy-momentum tensor as 
the sum 



D-l 



{Ti)p,q = (if )FRW + C^p.gj (T^l = ^ ^3+^ T i)j,D-j-, 



(53) 



where (T* )frw is the part corresponding to the uncompactified FRW spacetime and (Xf ) p *q is 
induced by the the nontrivial topology. The first term is well investigated in the literature (see 
[TJ [2] and references therein) and in the following we will concentrate on the topological part. 
Note that this part has the general structure (Zf)p*g = t~ D ~ 1 5^ /W(L p +i/r/, . . . ,Ld/t]), where 
the form of the functions f^\xi, . . . ,x q ) directly follows from formulae (I47h and (|49h . Recall 
that the ratios in the arguments of these functions may also be written in the form (|35p . 

The expressions for the VEV of the energy-momentum tensor are further simplified for a 
conformally coupled scalar field. In this case, after some calculations we find (no summation 
over i): 



A P+ i(T ? 



2ft 



-D-l 



(2vr)P/2+ 1 K x ^ E 



E 



006(27^0 {nL fcn _ l}) (54) 
(nL p+ i)P+ 2 ( c )p V p+ 9 ' K ' 
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where 

f{c)p( X ) = fp/2+l( x )i i = 0,l, 

= -(P + l)f P /2 + i(x)-x 2 f p/2 (x), (55) 

/(c)p( X ) = ( nL p+i k i) 2 f P /2(x), i = p+2,...,D. 

This result is also obtained from the corresponding formulae in the toroidally compactified 
Minkowski spacetime by using the standard relation for the VEVs in conformally related prob- 
lems. Formula (I54D for a conformally coupled field holds for any conformally flat bulk with 
general scale factor O(r). As wee see, for a conformally coupled field the stresses along the 
uncompactified dimensions are equal to the vacuum energy density. 

In the special case of power law expansion with c = 2/(D + 2) and for general £ we again 
have v = 1/2. After the integration over y, in this case we find 

,-D-l 



■ _ 2Q \- \- cos(27rna p+1 ) 

2 



(56) 



where fo = D and /j = 1 for i = 1,2, ... ,D. Recall that this special case corresponds to the 
radiation driven expansion and the corresponding Ricci scalar is zero. 

The formulae given above for the VEV of the energy-momentum tensor in the general case 
of the curvature coupling are simplified in the asymptotic regions of the ratio L p+ i/rj. For small 
values of this ratio the arguments of the modified Bessel functions in the expressions for the 
VEV of the energy-momentum tensor are large. In the leading order we have F(°)(z) « -z 2 ^ 1 
and to this order the first term in the square brackets in the expression (I49p for the vacuum 
stresses does no contribute. The integral over y is explicitly evaluated and we find 



(57) 

As we see, in the leading order the VEV coincides with the corresponding expression for a 
conformally coupled field. This limit corresponds to late stages of the cosmological expansion 
(t — > oo) for < c < 1 and to early stages (t — > 0) for c > 1. In terms of the comoving time 
coordinate, the topological part behaves as A p _|_i(T?) Pj(? oc £ _ (- D + 1 ) c . By using the relation (f35|) . 
we see that the asymptotic expression ([57|) corresponds to small proper length of the compact 
dimension compared to the Hubble length. 

For large values of the ratio L p+ i/r] and in the case of positive v for the topological part in 
the VEV of the energy-momentum tensor one has (no summation over i) 

kp+x{Tt)p, q « j^A p+1 (v 2 ) P , q , L p+1 /r) » 1, (58) 
where the asymptotic expression for the field squared is given by Eq. (|38|) and we have defined 
F<») = »0-„)» + 2nec ^ + <^p£, 

F«> - 2 ( t - i) 0, - „f - 2^ - lo ^-fS^^ , (59) 



12 



with I = 1, ... ,D. As it is seen, the leading terms in the vacuum stresses are isotropic. In 
terms of the comoving time coordinate, we have the behaviour A p _|_i(T?} Pj ,j oc t 2 ( c - 1 ) l '- cD — 1 m 

Note that this holds for the total topological part (T?) p *q as well. The limit under consideration 
corresponds to late times of the cosmological expansion (t — > oo) for c > 1 and to early times 
(t — > 0) for < c < 1, or alternatively, to large comoving length of the compact dimension in 
units of the Hubble length: L^^/th 3> 1. In the case v = the asymptotics of the topological 
part in the VEV of the energy- momentum tensor are still given by relation (|58p , where now 

/ = 0. (60) 



f0 , _ jD{D - 1) {l) = J_ 
' (1 - 1/c) 2 ' 2 



D - 1 

1 + 



(1 - 1/cf 



Note that the coefficient in (|58p does not depend on p and, hence, similar relation takes place 
between the total topological parts in the VEVs of the field squared and the energy-momentum 
tensor (given by formulae (|34p and (|53p ). In the limit under consideration the vacuum stresses 
are isotropic and the equation of state for the topological parts in the vacuum energy density 
and pressures is of the barotropic type (no summation over i): (T^pjj (F^ / F^)(Tq)p^, 
i = 1,2, . . . , D. In the case of a minimally coupled scalar field and for 1/D < c < 1 we have 
v = — b and, as it follows from (j59j) . A p+ i(T/) Pj(? ps — A p+ i(Tq)p j(? , i = 1, 2, . . . , D (no summation 
over i). This corresponds to the equation of state for stiff fluid with pressure equal to the energy 
density. 

For a minimally coupled scalar field in the case c £ (0, 1/D) U (1, oo) (in this case v = b) and 
for a conformally coupled field the functions FW in (|59p vanish. In these cases we should keep 
the next terms in the corresponding asymptotic expansions. For a conformally coupled field the 
behavior is described by the exact formula (|54p . In the case of a minimally coupled field and for 
c G (0, 1/D) U (1, oo), to the leading order we have (no summation over i): 

A p+1 (T > p , g \+l(Ti) p , q 2 b + 1 T(6)^0- 2 ^ cos(2vrna p+1 ) 



6-1 ~ b-2 (2vr)(P+ 3 )/ 2 y g _ 1 ^ {nL p+1 )P+ 3 - 2b 

x Yl /(p+3)/2-6(ra£p+ifcn,_i), * = (61) 

n^ieZ'- 1 

and for the stresses along the compact dimensions 

+1 ^ P+1-& A , t0 , , 2 fe + 1 r(fr)^- 2 ^ cos(27rna p+1 ) 
*P+lVp+l)p,i ~ b _ 1 *P+iVoh, g + ( 2vr )( P +3)/2y g _ 1 ( n L p+1 )P+2-26 

x ^ (nL p+ iA; nij _i) 2 /(p+i)/2-b(^p+i^n q _i), 
^ 2 fc +!r(6)^- 2 ^ cos(2vrna p+1 ) 



(2 ^ ) (p+3)/2y g _ 1 (nL p+1 )P+2- 



2/< 



x ^ {nL p+1 k i ff^ +1 y 2 -b{nLpj r ik aq _ 1 ),ip-\-2,...,D. (62) 

Hence, in this case the topological part behaves as A p+ i(r?) Pj9 oc t~ 2c in the limit t — ► oo for 
c > 1 and in the limit t — > for < c < 1/D. 

For small values of the ratio r]/L p+ i and imaginary v we use the asymptotic formula (|39j) . 
To the leading order this gives (no summation over I) 

2 (l-p)/2 Bn -D-l 

~ 7T ( P +3)/2y ^p ~ cos [ 2 IH M-^p+i/ 7 ?) ) Lp+i/rj^l, (63) 
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where the constants B\ and <pi are defined by the relation 

B^=2^T^\)F^ £ E ^^ /M^^VO. (64) 

ng-ieZ 1 ?- 1 n=l 

with F^, / = 0,1,..., D, from (JES). Hence, for imaginary f the behaviour of the VEVs is 
oscillatory. The corresponding amplitude behaves as t~ cD ~ l and we have damping oscillations. 

Summarizing the analysis given above, for a minimally coupled scalar field in the models 
with < c < 1, at late stages, t — > oo, the topological part in the VEV of the energy- momentum 
tensor behaves as t~( D+1 ^ c . At early stages, t — ► 0, we have ^Tf)p\ oc t~ 2c for < c < 1/D 
and {T} t )p\ oc t~ 2Dc for 1/D < c < 1. In the models with c > 1 and at late stages of the 
cosmological expansion the topological part decays like t~ 2c . At early stages and for c > 1 one 
has {Tl) { p\ oc t~( D+1 ) c . 



5 Special case and numerical examples 

As an application of the general formulae given above let us consider the special case of topology 
R ^ 1 x S 1 . For this case the expression for the topological part in the VEV of the field squared 
takes the form 



4(2vr)- D / 2 - 1 ^ cos(2irna D ) 



2x(t) 

' D - 1 ' 1 (Ut]) - 1 ^ (nL D /r]) D - 2 

/•oo 

x dyy[I- u (y) + I u (y)\ K u (y)f D/2 ^ l (nyL D /r ] ). (65) 
J o 

Note that in this formula Vt-q = t/\\ — c\. For a conformally coupled field this formula is reduced 
to 



/ 2 N (t) T((£> _ i)/ 2 ) cos(2irna D ) 

& )d-i,i ~ 27T (D +m{QLD)D -i n D-i > e - (66) 

n=l 



Formula (|66p describes the asymptotic behavior of the VEV in the general case of £ in early 
stages of the cosmological expansion (t — > 0) for c > 1 and at late stages (t — ► oo) for < c < 1. 
In the opposite limit of late times (t — > oo) for c > 1 and early times (t — > 0) for < c < 1, for 
positive values of u, the asymptotic behavior is given by the expression 

(v 2 ) {t) ~rfJ/2 t Q £M(W!?)___! y^ cos(27rna D ) 

V ' 71=1 

Note that in this expression QLd is the comoving length of the compact dimension. For imag- 
inary values of v we have the asymptotic described by the relation (|63p . In particular, in the 
models of power law inflation (c > 1) we have {f 2 )^_ l l oc t^~ D ^ c at early times. At late times 
the topological part behaves monotonically as j2(b—v)(i-c) j n ^ e case Q f positive v and oscillatory 
like i}~ cD cos[2|z^|(c — l)lni + const] for imaginary values of v. In the latter case the distance 
between the nearest zeros linearly increases with increasing t. 

Note that the ratio (cp 2 )^_ 1 i/(<P/i)_i 1 , with ((p 2 )^'^ 1 being the VEV for a conformally 
coupled field given by (IBB]) , is a function of Ld/t] alone and depends on the parameters £ 
and c through v. In figure [U we have plotted this ratio as a function of L/rj, with L = Ld 
being the length of the compact dimension, for untwisted D = 3 scalar field {ocd = 0) and 
for various values of the parameter u (numbers near the curves). In this special case one has 



14 



(v? 2 )£)_i i = l/[12(r2L£)) 2 ]. Note that the ratio L/rj is related to the comoving length of the 
compact dimension, measured in units of the Hubble length, by Eq. (135j) . Figured] clearly shows 
that the adiabatic approximation for the topological part is valid only for small values of the 
ratio L( c )/Vh. 




7 2 3 4 5 



Figure 1: The ratio (tp 2 )^_ l i/ { l f 2 )^^i i versus L/r) for different values of the parameter v 
(numbers near the curves) and for D = 3 scalar field with the periodicity condition along the 
compact dimension. 



In the special case under consideration, for the topological part in the VEV of the energy- 
momentum tensor we have the following representation (no summation over i) 



A(2tt)- d / 2 - 1 ^ cos (27rna D ) 



(Or?) 



D- 



— E 



; (nLn/r,)^ J 



dy y 



3-2b 



FW(y)fn/ 2 -i(nyL/r,) - ^^ f^nyLnh) 



(68) 



where F {0 \y) is defined by (@S]), f$_ x (x) = 0, and F®{y) = F(y) for i = 1, . . . , D. The 

functions f^^x) for i = 1, . . . , D are given by expressions (|50p . For a conformally coupled field 
we have (no summation over % = 0, 1, . . . , D — 1) 



1 l T D\(t) 



r((D + i)/2) 



^(D+l)/2(^ L 



D 



n=l 



cos(27rno;D) 



n 



D+l 



(69) 



Formula (|69|) describes the asymptotic behavior of the VEV for the general case of the curvature 
coupling parameter in the limit In the opposite limit the corresponding asymptotic 

formulae directly follow from the expressions (|58p and (j63H for real and imaginary u respectively. 
The formulae are simplified in the case c = 2/(D + 1) for general curvature coupling parameter 
(no summation over i): 



ctf>8 } -i,i 



T((D - l)/2) ^ CO s(27rnaz)) 



tt( d + 1 )/ 2 (OL d )^ 



n=l 



D+l 



D - 1 



£> - 1 



(70) 
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where gi = 1 for i = 0, 1, . . . , D — 1, gjj = —D and /j is denned after formula ([561) , 

In order to see the relative contributions of separate terms in the VEV of the energy- 
momentum tensor, in figure [2] we have plotted the ratio (T^)^_ 1 1 /(Tq )frw i n the radiation 
driven model (c = 2/(D + 1)) for D = 3 minimally coupled untwisted (left panel) and twisted 
(right panel) scalar fields versus rj/L with L = Lq being the length of the compact dimension. 
Here, (Tq)frw = l/[4807r 2 (2t) 4 ] is the vacuum energy density for a massless minimally coupled 
scalar field in the corresponding FRW model with trivial topology. Note that for the special case 
under consideration one has r\ = (2/a)t 1 ^ 2 . In the same model and for a conformally coupled 
scalar fields we have (no summation over i = 1,2) 

(T»)«'»=(T-)«'»=-i(T|)J>=0(|) 4 , (71) 

where C = —1/90 and C = 7/720 for untwisted and twisted fields respectively. 




0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 



Figure 2: The topological parts in the VEV of the energy-momentum tensor for the case of 
radiation driven model with spatial topology R 2 xS 1 as functions of the ratio rj/L. The left/right 
panel corresponds to untwisted/twisted scalar fields. 

Now let us estimate the relative contributions of different sources in the total energy density 
for the D = 3 radiation dominated model. Introducing the Planck mass Ml = 1/G and the 
Planck time tp, for the energy density being the source of the background metric one has 
p ~ Mp/(t/tp) 2 . The average vacuum energy density in FRW model with trivial topology is 
estimated as (Tq)frw ~ Mp/(t/tp) 4 . For the topological part in the vacuum energy density in 
the case of the conformally coupled field we have (^o )conf ~ 

M*/(L( C >/I P ) 4 , where = a(t)L 
is the comoving length of the compact dimension and lp is the Planck length. For a minimally 
coupled field one has (To)^ in ~ (TQ)^ nf [l + const- (L^ /t) 2 ]. From these relations, for the ratios 
of the separate parts in the vacuum energy density to the source term we find 

(T °) F rw/p ~ (t P /t) 2 , <T °)« f /p ~ (iprn/L^ 2 ) 2 , (72) 

where rn is the Hubble length. From these relations we conclude that the term (Tq)frw is 
comparable to the source term p only near the Planck time. The topological part dominates the 
part (Tq)fkw when < rn and becomes of the order of the source term for L^ 2 ~ Iprn. In 
the latter case the backreaction of quantum topological effects should be taken into account. 
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6 Conclusion 



In the present paper we have investigated one-loop quantum effects for a scalar field with general 
curvature coupling, induced by toroidal compactification of spatial dimensions in spatially flat 
FRW cosmological models with power law scale factor. We treat gravity as a given classical 
background field and do not consider the backreaction of the quantum effects on the metric. 
General boundary conditions with arbitrary phases are considered along compact dimensions. 
As special cases they include periodicity and antiperiodicity conditions corresponding to un- 
twisted and twisted fields. The boundary conditions imposed on possible field configurations 
change the spectrum of vacuum fluctuations and lead to the Casimir-type contributions in the 
VEVs of physical observables. Among the most important characteristics of the vacuum state 
are the expectation values of the field squared and the energy-momentum tensor. Though the 
corresponding operators are local, due to the global nature of the vacuum these VEVs carry an 
important information on the global structure of the background spacetime. As the first step 
we evaluate the positive-frequency Wightman function. By using the Abel-Plana summation 
formula, we have derived recurrence relation connecting the Wightman functions for the topolo- 
gies KP x (S 11 ) 9 and R p+l x (S 1 )^ 1 . The repeating application of this relation allows to present 
the Wightman function as the sum of the function for topologically trivial FRW model and the 
topological part. The latter is finite in the coincidence limit and in this way the renormalization 
of the VEVs for the field squared and the energy-momentum tensor is reduced to that for the 
FRW universe with trivial topology. The latter problem is well investigated in the literature 
and we concentrate on the topological parts. 

The topological parts are given by formulae ([33]) and ([33]) for the field squared and by 
formulae ()4"T]) . ([39]) . ([53]) for the energy density and the stresses. In the case of a conformally 
coupled scalar field the problem is conformally related to the corresponding problem in flat 
spacetime with topology RP x (S" 1 ) 9 and the general formulae are simplified to (I36p and (I54p 
for the field squared and the energy-momentum tensor respectively. These formulae describe 
the asymptotic behavior of the VEVs in the general case of the curvature coupling in the limit 
Lj/ij C 1. In this limit the comoving lengths of the compact dimensions measured in units of 
the Hubble length are small and it corresponds to early stages of the cosmological expansion for 
< c < 1 and to late stages of the cosmological expansion when c > 1. The topological parts 
in the VEVs behave like {tp 2 )p\ oc t c ( 1_D ) and (T^)p/ q oc t~ c ( 1+D ) and the stresses along the 
uncompactified dimensions are equal to the vacuum energy density. 

In the opposite limit Li/t] 3> 1, corresponding to large comoving lengths of compact di- 
mensions compared to the Hubble length, the behaviour of the VEVs is qualitatively different 
for real and imaginary values of the parameter v defined by Eq. (|15|) . For real values, the 
asymptotics are given by formulae ([38]) and ([58]) and the topological parts in the VEVs be- 
have as {if 2 ) { pl oc t2(c-i)i/-cD+i and (2*)JJ) oc t 2( - c - 1 >- cD ~ 1 . In this limit the vacuum stresses 
are isotropic and the equation of state for the topological parts in the vacuum energy den- 
sity and pressures is of the barotropic type. For a minimally coupled scalar field in the case 
c G (0, 1/D) U (l,oo) and for a conformally coupled field, the leading terms in the asymptotic 
expansions for the VEV of the energy-momentum tensor vanish. The corresponding asymp- 
totics are given by (Tfypjq oc £~ c ( 1+£) ) for conformally coupled field and by (2j*)p*q oc i _2c for a 
minimally coupled field in the model with c S (0, 1/D) U (1, oo) (see formulae ([6T]) and ([62]) ). In 
these special cases the vacuum stresses are anisotropic though the equation of state remains of 
the barotropic type with the coefficients depending on the lengths of compact dimensions. In 
the limit Li/r/ S> 1 and for imaginary values of the parameter v the asymptotic behaviour of 
the topological parts is described by the formulae ([40]) and ([63]) and this behaviour is oscillatory 
with the oscillating factor in the form cos[2|z/|(c — l)ln(i/io) + ip]- For the field squared the 
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amplitude of the oscillations behaves as t . For the energy-momentum tensor the amplitude 
of the oscillations decreases like t~ cD ~ l and the oscillations are damping for all values of the 
parameter c. 

The general results for the VEVs of the field squared and the energy-momentum tensor are 
specified in section [5] for the case of a model with a single compact dimension. By numerical 
examples for this case we have illustrated that the adiabatic approximation for nonconformally 
coupled fields is valid only in the limit when the comoving length of the compact dimension is 
very short compared to the Hubble length. We have also estimated the relative contributions of 
various sources in the total energy density. In particular, in the radiation driven model, for the 
comoving lengths < lprn the topological part of the energy density becomes of the order of 
the source term and the backreaction of quantum effects should be taken into account. 

In this paper we have considered the case of a massless field when the equation (jllH for the 
mode functions is exactly integrable in the class of special functions. By using this equation, 
we can estimate the role of effects due to the nonzero mass of the field. The ratio of the first 
and second terms in figure braces of (jlip is of the order m 2 t 2 . Hence, for mt <C 1 the effects 
related to the nonzero mass may be disregarded and the results described in the present paper 
are applicable for massive fields as well. 
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